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Abstract. Recent X-ray observations of the rich galaxy 
cluster CL 0939+4713 (or Abell 851) at z = 0.41 with the 
ROSAT/HRI as well as with the ASCA/GIS and ASCA- 
/SIS instruments are presented and analysed. With the 
high resolution imaging data (ROSAT/HRI) we confirm 
and extend the earlier impression that the cluster has 
a lot of substructure. Two maxima of the cluster emis- 
sion are obvious in the images, each of them shows even 
some internal structure. One of the subclusters can be 
nicely modeled with an elliptical model. For the total lu- 
minosity of the cluster in the ROSAT band we obtain 
Lx,o.i~2AkeV = 6.4j^3 x 10 44 erg/s, for the bolomet- 
ric luminosity Lx.boi = 1.6jZoli x 10 45 erg/s. We perform 
spectral fits for the two ASCA instruments and for the 
ROSAT/PSPC simultaneously. The most reliable num- 
bers for the temperature and metallicity for the intra- 
cluster gas are T = 7.6 j^g keV and m — 0.22^q;22 to ©' 
respectively. We find a relatively small total mass, a small 
gas mass ratio and a small iron mass to light ratio. These 
numbers together with relatively low luminosity for such 
an optically rich cluster and the pronounced substructure 
confirm that CL 0939+4713 is a young cluster still far 
away from a virialised state. 

In the same ROSAT/HRI image the X-ray emission of 
a background quasar (zq w 2) can clearly be identified. 
With a luminosity of Lxfi.x— lAkeV ~ 

1.4 x 10 45 erg/s 

(which is necessarily affected by gravitational lensing) it 
belongs to the X-ray brightest quasars. A striking differ- 
ence between the ROSAT/PSPC and the ROSAT/HRI 
observations (taken almost five years apart) is the faint- 
ing of one X-ray source (RXJ0943. 0+4701) by at least a 
factor of ten. We try to identify this source from deep op- 
tical data. Our best candidate is a blue compact object, 
possibly an AGN. 
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1. Introduction 

One of the foremost goals of modern cosmology is the un- 
derstanding of the formation and evolution of structure. 
Clusters of galaxies provide an excellent probe because 
they represent the largest gravitationally bound objects 
in the universe and can be identified over a large range of 
redshifts. Studying the X-ray properties of distant clusters 
and comparing them with the properties of nearby clus- 
ters gives important insight into the evolution of galaxy 
clusters and provides a strong test for cosmological mod- 
els. 

The galaxy cluster presented in this paper is the rich 
cluster CL 0939+4713 at redshift z = 0.41. It is very well 
studied optically (see e.g. Dressier & Gunn 1983; Dressier 
& Gunn 1992; Dressier et al. 1993; Belloni et al. 1995; Bel- 
loni & Roser 1996). Almost 200 (spectroscopic and photo- 
metric) redshifts are known. The cluster is extremely rich, 
and it contains a high fraction of blue and E+A galax- 
ies. A background quasar at redshift two is seen through 
the cluster. CL 0939+4713 was subject of a weak lensing 
analysis (Seitz et al. 1996), and very recently Trager et al. 
(1997) showed that there are a few redshift four galaxies 
seen through the cluster. 

A ROSAT/PSPC observation of this cluster was analy- 
sed by Schindler & Wambsganss (1996). It was shown 
that the X-ray luminosity is relatively low for such a rich 
cluster. Meanwhile new X-ray data are available from the 
high resolution ROSAT/HRI detector as well as improved 
spectral information from the ASCA/GIS and SIS instru- 
ments. 

Here we present and analyse this new X-ray data. We 
investigate the X-ray morphology of the cluster (Sect. 2.1 ) 
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from the new ROSAT data and show that the cluster con- 



tains clear substructure. In Sect. 2.2 we present a spec 



tral analysis from both ASCA instruments and from the 
ROSAT/PSPC data. The X-ray luminosity is determined 
in Sect. [2.3|. In Sect. 2.4 we estimate the gas mass, the iron 



mass and the total mass. We look in some detail into the 
variable X-ray source RXJ0943. 0+4701 whose flux weak- 
ened by more than a factor of ten over the five years be- 
tween the two ROSAT observations (Sect. |2^). The at- 
tempt to identify an optical cou nterpart for RXJ0943.0- 
+4701 is described in Sect. 2J3. Finally we identify and 
measu re t he X-ray properties of the quasar in the field 
(Sect. |2.7|) . In the final Chap. || we discuss and summarize 
the X-ray and optical properties of this unusual galaxy 
cluster. 

Throughout this paper we use Hq = 50 km/s/Mpc. 
2. Observations 

For the analysis we use optical and X-ray observations - 
taken with ROSAT (Triimper 1983) and ASCA (Tanaka 
et al. 1994) - obtained within a period of five years. The 
X-ray observations consist of three data sets taken with 
different instruments: 

- 14.4 ksec with ROSAT/PSPC on November 15 th to 
20 th , 1991 

- 39 kscc with ASCA on April 12 th and 13 th , 1995 

- 45.6 ksec with ROSAT/HRI on October 16 th to 20 t/l , 
1996 

While the ROSAT observations give us information on 
the spatial distribution of the X-ray emission, the ASCA 
data - and to some degree also the ROSAT/PSPC data - 
provide spectral resolution. The ASCA analysis uses both 
types of detectors - the Solid-state Imaging Spectrometers 
(SIS) and the Gas Imaging Spectrometers (GIS). The SISs 
were operated in 2CCD Faint mode, while the GISs were 
operated in PH mod e (see http: / / asca.gsfc.nasa.gov/docs 



/asca/abc/abc.htm] ). The ASCA data are selected ac- 



cording to the following criteria: elevation higher than 20 
degrees and bright Earth angle larger than 30 degrees. 

The optical observations have been all obtained in 
November 1992 at the prime focus of the 3.5 m telescope 
on Calar Alto (Spain) and are listed in Table 1 of Bel- 
loni et al. (1995). A series of narrow-band filters had been 
used (FWHM ~ 90-200 A) covering the range from 3800 
to 9200 A, and broad band B, R, I filters to obtain low res- 
olution spectra for all galaxies brighter than R=22.5 mag 
in a 7 arcmin x 7 arcmin cluster field. The filters were 
selected to examine as accurately as possible the spectral 
energy distribution (SED) of ellipticals and E+A galaxies 
(Dressier & Gunn 1983) at the clusters' redshift. There- 
fore, they allow us to estimate the redshift of these galax- 
ies with an accuracy a z = 0.010 (Thimm et al. 1994) and 
with an accuracy of about a z = 0.030 for cluster members 



of other morphological types and foreground/background 
galaxies. 

2.1. X-ray morphology of CL 0939+4713 

The ROSAT/PSPC data of cluster CL 0939+4713 were 
analysed and discussed in Schindler & Wambsganss 
(1996). Here we concentrate on the higher resolution and 
much longer HRI observation and compare the two data 
sets. 

The ROSAT/HRI image of CL 0939+4713 obtained in 
October 1996 is shown in Fig. 1. The cluster looks very 
amorphous. It it obvious that there is a lot of substructure 
in this cluster. The cluster emission consists mainly of two 
maxima labeled Ml and M2 (see Fig. 1.). Both of them 
seem not spherically symmetric but show some internal 
structure. The positions of these maxima are: 

- Ml at 9 h 42 m 58.2 s , 46°58'52" (J2000), with a maxi- 
mum surface brightness of 9.2 x 10~ 3 erg/arcmin 2 /s, 
and 

- M2 at 9 h 43 m 04.4 s , 46°59'48" (J2000) with a maximum 
surface brightness of 10.2 x 10 -3 erg/arcmin 2 /s. 

These values are obtained from an image smoothed with 
a Gaussian filter of a — lOarcsec. The central part of the 
ROSAT/HRI image looks very similar to the smoothed 
light distribution of all identified E/S0 galaxies (Seitz et 
al. 1996, Fig. 9b) as observed by Dressier et al. (1994). 

The high resolution of the HRI allows for the first time 
to resolve the emission of the quasar at z=2.055 (Dressier 
et al. 1993). It is located left of the left maximum and 
marked with a plus sign in Fig. 1. 

The most striking difference between this HRI ob- 
servation and the earlier PSPC image is the complete 
absence of the emission in the north which is marked 
as M3 (RXJ0943.0+4701) in the PSPC observation (see 
Schindler & Wambsganss 1996). This will be investigated 
in detail in Section ^1]. The bright point source (PI) in 
the left part of Fig. 1 was detected already with the PSPC. 
It is very likely not related to the cluster. 

As the cluster is clearly not spherically symmetric we 
do not attempt to fit a single (3 model (Cavaliere & Fusco- 
Femiano 1976; Jones & Forman 1984) 



S(r) = S 1 + 




-3/3+1/2 



(1) 



to the cluster as a whole, but to the two subclusters sep- 
arately (S(r) is the surface brightness at distance r, So is 
the central surface brightness, r c is the core radius, and 
(3 is the slope parameter). The fits are centred on the re- 
spective maximum, a sector of 90° in the direction of the 
other subcluster and the quasar emission being excluded. 
The profiles and the fits are shown in Fig. 2, the fit pa- 
rameters are listed in Table 1. Both subclusters show very 
small values for (3 and for the core radius. 
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09 h 43 m 20 s 09 h 43 m 00 s 09 h 42 m 4( 

Fig. 1. ROSAT/HRI image of the cluster CL 0939+4713. The image is smoothed with a Gaussian filter of a = 
lOarcseconds. The cluster shows very pronounced substructure. The two main maxima are marked with Ml and M2. 
The extension left of the second maximum is caused by a quasar at z = 2.055. Its exact position is indicated by the 
plus sign. The bright point source (PI) is very likely not related to the cluster. 



For a quantitative analysis of the substructure we ap- 
ply the method by Neumann & Bohringcr (1997): we fit an 
elliptical (3 model to the HRI data. Since the cluster has 
a bimodal shape an elliptical fit cannot converge. There- 
fore we have to exclude one of the maxima for fit. Both 
maxima have about the same number of counts within 
a radius of 40 arcseconds, but Ml is less peaked and its 
shape corresponds less to an elliptical (3 profile. There- 
fore we exclude Ml and centre the fit on M2. The quasar 
emission is excluded for the fit as well. The resulting fit pa- 
rameters are given in Table 1. They correspond well with 
the parameters of the spherically symmetrical fit, hence 



we can be sure that the elliptical fit converged reasonably 
well. Figure 3a shows the elliptical model. 

Subsequently, we use the elliptical model and subtract 
it from the original image. The residuals are shown in 
Fig. 3b. The emission of and around M2 has almost dis- 
appeared which indicates that the elliptical fit is a good 
model. 

In this residual map the maximum Ml is clearly visible, 
it has a significance of more than 4cr. It is quite obvious 
from the morphology that Ml is not relaxed but has some 
internal structure. There is also some diffuse emission sur- 
rounding Ml, but with lower significance. The emission 
of the quasar is very clearly visible. The other emission 
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r (arcmin) r (arcmin) 

Fig. 2. Surface brightness profiles around the maxima Ml and M2 with a spherically symmetric /3-modcl fit. The fit 
parameters are given in Table 1. 



Type of fit 


centre 


So in cts/s/arcmin 2 


r c in arcsec 


r c in kpc 


P 


PA (N over E) 


spherical fit 


Ml 


5.8 x 1(T 3 


15 


95 


0.31 




spherical fit 


M2 


7.9 x 1(T 3 


11 


67 


0.36 




elliptical fit 


M2 


8.3 x 10~ 3 


12 (major axis) 


79(major axis) 


0.36 


133° 








9(minor axis) 


56(minor axis) 







Table 1. Fit parameters for the spherically symmetric and elliptical (3 model fits to the two subcluster: So is the 
central surface brightness, r c is the core radius (for the elliptical fit we give two numbers for the major and the minor 
axes), (3 is the slope parameter, and PA is the position angle (only for the elliptical fit). 



around the cluster has only significances between 2 and 
3cr, i.e. it is not clear whether each of these features is 
indicating more substructure is real or whether some of 
them are just statistical fluctuations. 

2.2. Spectral analysis 

In order to complement the morphological analysis of CL 
0939+4713 we perform a spectral analysis as well. For that 
we use data obtained by both ASCA/GISs (G2 and G3) 
and ASCA/SISs (SO and SI) detectors and the ROSAT- 
/PSPC. 

For composing the ASCA spectra the photons are ex- 
tracted from a circle centred on 9 ,l 42 m 55.5 s , 46°59'47" 
with a radius of 3.5 arcminutes and 4.5 arcminutes for 



the SIS and GIS, respectively. It is necessary to extract 
the data from different radii because of the different point 
spread functions of these instruments. 

The GIS spectrum covers an energy range of 0.7-10 
keV and is rcbinncd to contain at least 50 photons per 
bin. The background is taken from blank sky fields where 
the contaminating sources are masked (Ikebe 1995). 

The SIS data are difficult to analyse as they were 
taken with a non-standard energy discrimination thresh- 
old. While normally a lower energy discrimination thresh- 
old of 0.41 keV is used, these SIS data were taken with 
a threshold of 0.55 keV. The correction of the detector 
response seems to be inaccurate for this non-standard 
threshold, e.g. if we use the SIS energy range down 
to 0.7 keV, we find very high tie values not consis- 
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tent with the PSPC results. Therefore we have to re- 
strict the SIS energy range for the spectral fits. In Ta- 
ble 2 we list the fit results for two different SIS energy 
ranges 3-6 keV and 1-10 keV. Obviously, these are both 
valid energy ranges because they give very similar re- 
sults. The SIS data are rebinned to contain at least 20 
photons per bin and the standard Godard background 



responding 90% confidence level errors: T = 7.&t\% keV, 



see http:/ /asca. gsfc.nasa.gov/docs/asca/abc/abc. html) 



is used. In the ASCA images a possible contribution of 
RXJ0943. 0+4701 is not resolved, therefore this emission 
can affect our results. But as RXJ0943. 0+4701 has a softer 
spectrum than the cluster (see Sect. 2.5), its influence 
is expected to be very small on the ASCA spectra. Ex- 
trapolating the flux of RXJ0943. 0+4701 as measured in 
the ROSAT/PSPC observation to an energy band of 1- 
10 keV - assuming a power law spectrum with a mean 
photon index of 2.5 (see Sect. 2.5) - we find that the flux 
-Fkxj0943(l-10keV) is less than 1% of the cluster flux. 

The source spectrum of the ROSAT/PSPC is ex- 
tracted from a circle centred between Ml and M2 
(9 h 43 m 01.6 s , 46°59'40") with a radius of 2.3 arcmin- 
utes (the point spread function is even smaller than 
the one of the ASC A/SIS). We exclude the emission of 
RXJ0943. 0+4701 (9 h 42 m 56.8 s , 47°00'47" with a radius 
of 50 arcseconds) since it is most likely unrelated to the 
cluster. The background is taken from different parts of 
the same pointing. The spectrum is rebinned to contain 
at least 16 and 25 photons, respectively. Fits with dif- 
ferent background realizations and different binning give 
very similar results. 

We fit the spectra of all instruments simultaneously by 
folding for each spectrum the Raymond-Smith model with 
the response for the corresponding instrument (see Fig. 4). 
This procedure has the advantage of enabling us to con- 
strain different fit parameters at the same time. While the 
PSPC data can determine the temperature only with very 
large uncertainty and cannot constrain the metallicity at 
all, they can determine the hydrogen column density quite 
reliably, because the nn affects mainly the soft part of 
the spectrum. The CIS data are ideal for determining the 
temperature, because the spectrum extends to the highest 
energies (see Fig. 4). The SIS data are the ones that give 
good constraints on the metallicity from the Fe K emis- 
sion line, as the energy resolution of the SIS is the best 
of all used detectors. Because all the fit parameters are 
connected, it is the best to fit the spectra of all detectors 
simultaneously. 

The results of the spectral analysis are listed in Table 
2. We show different simultaneous fit results, for two SIS 
energy ranges, each one with the redshift as free fit pa- 
rameter as well as the redshift fixed to the optical value 
of z — 0.41. The results of all fits agree well within the 
errors. The hydrogen column density is consistent with 
the Galactic value of 1.27xl0 20 cm -2 (Dickey & Lock- 
man 1990) within the uncertainty. In the following we use 
the mean values from all four simultaneous fits and cor- 



m = 0.221^7710 and n H = 0.9+^ x 10 20 c m -+ 

We find small discrepancies in the fluxes of the differ- 
ent detectors: the flux determined by GIS is about 15% 
higher with respect to the flux determined by SIS and 
about 25% higher with respect to the PSPC flux, which 
can be explained only partially by the different extraction 
radii and variability of the source RXJ0943. 0+4701 (see 
Sect. |J). 

For comparison with the temperature found previously 
in a ROSAT/PSPC analysis by Schindler & Wambsganss 
(1996) we show in Table 2 also the fit of the PSPC data 
only. The temperature of T= 5.7^^° 7 keV is somewhat 
lower than in the combined fits, but the upper limit is 
unconstrained. The temperature of Schindler & Wambs- 
ganss (1996) T = 2.9lg j keV obtained without exclud- 
ing RXJ0943. 0+4701 was much lower. The explanation is 
that the ROSAT/PSPC spectrum was contaminated by 
the soft contribution of the variable source RXJ0943.0- 
+4701 and therefore the best fit yielded a lower tempera- 
ture. 



2.3. Luminosity 

The X-ray emission of the galaxy cluster can be traced 
out to 2.5 arcminutes (almost 1 Mpc) in the HRI observa- 
tion. Within this radius we find 1060 source counts after 
background subtraction (the point source PI lies outside 
this circle; the emission of the quasar is excluded.) This 
corresponds to a count rate of 0.023 HRI counts/s. 

Using the temperature T — 7.6^ 2 'g keV and the hy- 
drogen column density nn = 0.9±q4 x 10 20 cm~ 2 of the 
spectral fits we derive a luminosity in the ROSAT band 
(0.1-2.4 keV) of 6.4±° j x 10 44 erg/s, corresponding to a 
bolometric luminosity of 1.6io'g x 10 45 erg/s. 

2.4- Gas mass, total mass and iron mass 

A mass determination from the X-ray data for this clus- 
ter is very difficult, because one has to assume spherical 
symmetry, which is certainly a bad approximation for CL 
0939+4713. Just for giving a rough estimate of the mass 
we restrict the mass analysis to the central parts of the 
two subclusters assuming that each subcluster is spheri- 
cally symmetric. The outer radius of 41 arcseconds (265 
kpc) is chosen as half of the distance between the subclus- 
ters. For the analysis we u se t he spherically symmetric 
/3-models described in Sect. 2.1 and a constant tempera- 
ture of 7.6 keV for both subclusters. 

Figure 5 shows the integrated gas and total masses. 
At a radius of 265 kpc the integrated gas masses in the 
subclusters Ml and M2 are very similar: MMi,gas(?" < 
265kpc) = 4.1 x 10 12 M Q and M M2>gas (r < 265kpc) = 
4.2 x 10 12 M©, respectively. Assuming a constant tem- 
perature and hydrostatic equilibrium we can also deter- 
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Instruments 


SIS energy range in keV 


k B T in keV 


m/rriQ 


n H in 10 20 cm -2 


z 


d.o.f. 


X 2 /d.o.f. 


GIS/SIS/PSPC 
GIS/SIS/PSPC 
GIS/SIS/PSPC 
GIS/SIS/PSPC 
PSPC only 


3-6 
3-6 
1 - 10 
1 - 10 


7.4± 2 ' 4 

8.i+i 3 7 

7.2+™ 
7 6+ 1 ' 6 

c 7+00 
O.I _ 2 .7 


C\ OO+0.23 

u -^°-0.20 

20+ 21 

u - zu -0.20 

2^ +0 - 21 

1Q+ - 20 

0.23(fixed) 


0.9±°oi 

0.8±° o j 

n Q+ - 5 
u.y_ .4 

n o+ - 5 

u - a -0.4 

1 o +0 - 6 


n oo+0. 07 
0.38_ 04 

0.41(fixed) 

n OO+0.05 

U-^-o.os 
0.41(fixed) 

0.41(fixed) 


81 
82 
139 
140 
18 


0.97 
0.97 
0.92 
0.92 
0.86 



Table 2. Results of the simultaneous spectral fits of the data from both ASCA/SIS detectors, from both ASCA/GIS 
detectors and from the ROSAT/PSPC. A fit of the ROSAT/PSPC data alone is shown as well. The columns show 
the temperature fceT in keV, the metallicity m in solar units, the hydrogen column density uh, the X-ray determined 
redshift z, the degrees of freedom of the fit and the reduced x 2 ■ The errors are 90% confidence level. For the uncertainty 
in redshift a systematic error of 1% in the ASCA energy calibration is included. 



mine the total masses of the subclusters: MMi.tot(?* < 
265kpc) = 6.0±l 3 6 x 10 13 M Q and M M2 ,tot(r < 265kpc) = 
7.4I2 Q x 10 13 M Q , respectively. The errors come mainly 
from the uncertainty in the temperature (see Sect. |2.2| ), 
which are t.22% • ^ n additional uncertainty is introduced 
by the internal substructure of the subclusters which tends 
to cause an underestimation of the mass (Schindler 1996). 
The errors on the mass include another uncertainty of 15% 
error coming from this effect, from possible temperature 
gradients and from possible temperature differences be- 
tween the subclusters. 

With these numbers we find gas mass fractions of (7± 
2)% and (6 ± 2)% for Ml and M2, respectively. These 
fractions are very small compared to other clusters where 
one finds usually gas mass fractions between 10% and 30% 
(Bohringer 1995). 

For a comparison with the total mass estimate from 
gravitational lcnsing we determine the mass in the region 
were the HST field and our two circles around the subclus- 
ter centres overlap. The HST field is much smaller than the 
HRI field so that only part of the two circles are covered 
by the HST field. As Ml and M2 are close to the border of 
the HST field we use only the sectors of the circles which 
have overlap with the HST field: angles +128° - +294° 
(N over E) centred on Ml and angles —90° — 1-90° cen- 
tred on M2. The total overlapping area is only about 1.4 
arcmin 2 . For the comparison we have to determine the 
surface mass density i.e. projecting the mass onto a plane. 
For doing this we have to extrapolate and project the (3 
model. As the outer integration radius in the line of sight 
we choose 1 Mpc, which is a compromise between extrapo- 
lating too far into non-relaxed regions and summing most 
of the mass. Integrating out to a radius of 1.5 Mpc would 
increase the mass only be 3.5% which is negligible consid- 
ering the error. With this procedure we find a projected 
mass of l.O^o 3 x 10 14 M Q in the overlapping area. The er- 
ror is again the same percentage as used for MMi.tot and 

-^M2,tot- 



The mass from gravitational lensing as calculated 
by Seitz et al. (1996) in the same overlapping area is 
2.3 x 10 14 M Q . For this number a mean redshift of the 
background galaxies of z = 1.0 is assumed. If one would 
use a different redshift z — 0.8 and z — 1.5, the mass 
would be 2.8 x 10 14 M Q and 1.8 x 10 14 M Q , respectively. 
As this redshift is the largest uncertainty in the lensing 
mass estimate, it remains a discrepancy of a factor 2-3 
between the X-ray and the lensing mass. 

For determining the iron mass Mp e in the intra- 
cluster gas we use the metallicity from the spectral fits 
m = O.22lg22"i0- With this we calculate the ratio of 
M Fe /Mg as = 3.2tlj> x 10~ 4 . For the determination of 
the gas mass we use the two spheres of 265 kpc radius 
around the two subclusters. The iron mass in this vol- 
ume is Mpe — 2.7 x 10 9 Af Q . Within the same radii in the 
optical field are 26 galaxies. Estimating the total blue lu- 
minosity including an extrapolation to faint galaxies we 
find about Lb = 6 x 10 11 L q . This results in an iron mass 
to light ratio M Fe /L B «5x lO _3 A/ /L Q . 

For comparison we do the same also in a larger sphere 
of 584 kpc around Ml (the maximum radius completely 
covered by the optical field) by extrapolating the /^-profile. 
The subcluster around M2 is neglected in this case. For 
this case we find Mp e — 7.7 x 10 9 M Q , 75 galaxies with 
L B = 15xl0 n L©, and again M Fe /L B sa 5xlO~ 3 M /L . 
As the two numbers for the iron mass to light ratio are so 
similar, we can be sure that we do not make an error by 
using a restricted volume and by not taking into account 
the gas in the line-of-sight outside the circles. 

The iron mass to light ratio Mf £ /Lb ~ 5 x 
10~ 3 Mq/Lq is smaller than the ones typical for rich clus- 
ters of 0.01-0.02 (Renzini et al. 1993), although one has to 
keep in mind the errors only coming from the metallicity 
are already ±100%. 
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2.5. The variable X-ray source RXJ0943. 0+4701 

A comparison of the PSPC and the HRI observations 
shows that the source RXJ0943. 0+4701 is not present in 
the HRI image (Fig. 6). For a determination of the flux 
of RXJ0943.0+4701 in the PSPC observation we simulate 
the PSPC point spread function, smooth it with a Gaus- 
sian of a — 20arcsec and subtract it from a PSPC image 
with the same smoothing. We vary the number of counts 
in the point spread function model until the level of the 
residuals is comparable to the diffuse cluster emission at 
this position. 

With this method we find that RXJ0943.0+4701 has 
about 50 counts (broad band) above the cluster emis- 
sion in the PSPC observation. This corresponds to about 
3.5 x 10~ 3 cts/s for the PSPC detector. With the above 
mentioned value of the Galactic nn and a power law with 
photon index 2.3 this count rate corresponds to a flux 
in the ROSAT band (0.1-2.4 keV) of / RO sat,Nov.91 = 
4 x 10~ 14 erg/cm 2 /s. 

In the HRI image appears very faint emission near the 
position of RXJ0943. 0+4701 with very low significance 
(only 1.3c) which could as well be a statistical fluctua- 
tion. If there is any real emission at this position it cannot 
contain more than 5 counts above the cluster plus back- 
ground emission; this corresponds to a source countrate 
< 1 x 10~ 4 HRI cts/s. With the same assumptions as above 
this count rate corresponds to a flux in the ROSAT band 
/ROSAT,Oct.96 < 5 x 10~ 15 erg/cm 2 /s. Given the smaller 
point spread function of the HRI which would give the 
(point) source a higher significance above the background 
and the diffuse cluster emission we conclude that the de- 
crease in luminosity of RXJ0943.0+4701 in the ROSAT 
band is at least a factor of 10. This decrease can be either 
due to a general decrease in luminosity or due to a shift 
of the emission out of the ROSAT energy range. 

The non-detection of RXJ0943.0+4701 in the HRI ob- 
servation cannot be caused by the fact that the HRI is 
slightly less sensitive in the soft band compared with the 
PSPC (ROSAT User's Handbook |http: / /www.xray.mpc" 
mpg.de/rosat/doc/ruh), because the hardness ratio at 
the position of RXJ0943.0+4701 (hard-soft)/ (hard+soft)- 
= —0.08, i.e. there are about the same number of counts 
in the soft and in the hard band. 

In the ROSAT All-Sky Survey (Voges et al. 1996) 
CL 0939+4713 was observed in November 1990 for two 
days. Unfortunately, the exposure time is too low to see 
any morphological details. The Survey countrate of the 
whole cluster region is within the errors in agreement with 
both the PSPC and the (converted) HRI count rate. A 
lightcurve of the Survey observation shows no indication 
of variability (Boiler & Voges, private communication). 

Unfortunately, RXJ0943. 0+4701 cannot be resolved 
with ASCA, i.e. one cannot distinguish the cluster emis- 
sion from emission coming from RXJ0943. 0+4701. There- 



fore, we have no information about the flux of RXJ0943.0- 
+4701 at the time of the ASCA observation. 

For examining the short-term variability of 
RXJ0943. 0+4701 we try to derive a lightcurve from 
the PSPC observation. The source was observed in eight 
intervals within 72 hours. The intervals have exposure 
times between 1300s and 2400s. Figure 7 shows this 
attempt of a lightcurve. It is consistent with a constant 
emission over 3 days. 

A spectral fit to the RXJ0943. 0+4701 region alone is 
very difficult because of the small number of photons. In 
an attempt to get a least a rough idea we fit the spectrum 
using two different kinds of background: 1) only detector 
background and 2) detector background plus cluster emis- 
sion. A fit with a Raymond & Smith (1977) model yields 
temperatures of T = l.lkeV and T = 0.3keV, respec- 
tively, for the two different background models. Although 
these results have very large errors they show that the 
RXJ0943. 0+4701 emission is considerably softer than the 
cluster emission (compare also Fig. 6). A fit with a power 
law to the same spectra yields photon indices of 2.2 and 
2.9, respectively (again with very large errors). 

2.6. Identification of RXJ0943. 0+4701 

Since the source RXJ0943. 0+4701 appears to be variable 
and hence very interesting, we try to find the optical coun- 
terpart of the X-ray emission at RXJ0943.0+4701. Unfor- 
tunately, the northern part of the cluster was never ob- 
served with HST. A ground based I-band image of the 
field around RXJ0943.0+4701 is shown in Fig. 8. The po- 
sition of RXJ0943. 0+4701 is not very well determined in 
the PSPC observation for several reasons. First, the max- 
ima in different bands are not exactly coincident (see Fig. 
6 and 8). Second, there is no other point source near the 
cluster with an obvious optical counterpart which could 
be used to correct for a possible pointing offset of the 
ROSAT telescope. Therefore, we use the HRI image for 
the correction. In this image the quasar at z=2.055 is used 
to determine the pointing offset of the HRI. The pointing 
offset is only 1 arcsec. After correcting for this offset we 
use the point source PI (see Fig.l) present in both, HRI 
and PSPC image, to correct for the PSPC offset, which is 
relatively large, 11 arcsec. 

The corrected positions of RXJ0943. 0+4701 are 
marked in the optical image (Fig. 8). In the region of these 
positions are 8 possible optical counterparts with R mag- 
nitudes brighter than 22 m . These objects are marked with 
numbers in Fig. 8. We measure the brightness of all these 
objects using fixed aperture photometry with an aperture 
of 3.5 arcsec on the deep cluster images taken with the 
12 filters (Belloni et al. 1995). The R magnitudes of all 
these counterpart candidates are listed in Table 3 along 
with their classification from the SEDs plus morphologi- 
cal analysis. Figure 9 shows the observed fluxes and the 
SEDs that we assign by eye as the best fit. Although this 
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Fig. 6. Morphology of CL 0939+4713 as seen in an HRI image (a) and in a PSPC image (b) taken 5 years earlier. 
The HRI image is smoothed with a Gaussian of a = 14arcsec so that the resolution is about equivalent to the PSPC 
observation, b) shows the PSPC image in broad band which is the corresponding band to the HRI image. The PSPC 
images are smoothed with a Gaussian of a = lOarcsec. In each of the panels the positions of the maxima of the 
PSPC broad band image are marked with crosses. The comparison shows that the northern maximum of the PSPC 
observation (RXJ0943. 0+4701) has disappeared in the HRI observation. The relative brightness of the maxima Ml 
and M2 seems to have changed, but this difference is within the expected statistical fluctuations, c) PSPC soft (0.1-0.5 
keV"l and d) PSPC hard CO. 5-2.0 keV") hand images a.re shown to demonstrate that RX.I0943. 0+4701 is softer than 
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is admittedly a quite qualitative classification, we can still 
gain valuable information from it. Using all this informa- 
tion we gather: of the possible optical counterparts two are 
stars, five are galaxies and one is a blue compact object 
(#130). 

The two stars, objects #166 and #115, have clearly 
very red spectra typical of cool stars. A convolution of 
K- and M-type stars spectra (Gunn & Stryker 1983) with 
our filter setup shows that both stars can be classified as 
K6-K8 or M0-M2 stars. A possible Ha emission would not 
be visible in the SEDs because the corresponding filter at 
6560 A has a FWHM of 170 A and therefore only a strong 
emission line (EW> 30 A) could be detected. But it is 
very rare that these stars have such is a strong emission 
line. 

For stars of this type an X-ray variability of a factor 
of ten is highly unlikely, unless the star is observed during 
a flare. But a flare can be excluded from the lightcurve, 
which shows approximately constant emission over 3 days 
(see Fig. 7). Furthermore, the ratio of X-ray flux to bolo- 
mctric optical flux would be too high for K or M stars. 
For the brighter star #166 the flux ratio 

/x,_RXJ0943.o+470i(0.1-2.4/ceV r )// optj# i 66 (6oO w 5xl0~ 2 , 

while a value of 10~ 3 is the upper limit for K or M stars (J. 
Schmitt, private communication). Therefore we can rule 
out the stars as origin of the X-ray emission of RXJ0943.0- 
+4701. 

It is also not likely that a cataclysmic variable causes 
the X-ray emission because a hardness ratio close to — 1 
(and not the observed —0.08) would be expected. 

None of the five galaxies shows the SED or the emission 
lines typical of an AGN. They look indeed like normal 
spiral galaxies, either in the foreground or possibly even 
cluster members. But we cannot exclude that one of the 
galaxies changed from an active state to a quiescent state 
within the year that passed between the PSPC and the 
optical observation. Even for an AGN a variability larger 
than a factor of 10 is quite exceptional, but an AGN would 
be a possible identification. 

As seen in the previous section, a spectral fit of the 
RXJ0943. 0+4701 region with a power law yields a photon 
index between 2.2 and 2.9 (with large errors because of the 
few photons and contamination from the cluster emission). 
These numbers are typical photon indices for AGNs (e.g. 
Boiler et al. 1996). 

Assuming that a cluster member galaxy (z = 0.41) 
or the galaxy #99 at a redshift of about z = 0.33 were 
the X-ray emitting AGN we calculate the luminosity in 
the ROSAT band. With a photon index of 2.4 and a 
Galactic uh we find a luminosity log[X(0.2keV)] = 26.5 
(corresponding to Lx(0.1-2.4keV)= 4.1 x 10 43 erg/s for 
z = 0.41) or log[L(0.2keV)] = 26.3 (equivalent to Lx(0.1- 
2.4keV)= 2.5 x 10 43 erg/s for z = 0.33). Both numbers 
are in the standard range for AGNs (Boiler et al. 1996). 



The values for Jx,rx jo943.o+47oi(0.1 - 2AkeV)/f opt 
are for all the galaxies in the typical range for AGNs (Mac- 
cacaro et al. 1988). So an AGN is a viable possibility for 
the source RXJ0943.0+4701. 

The most interesting candidate for the optical identi- 
fication of the strongly variable X-ray source RXJ0943.0- 
+4701, however, is the blue compact object #130. With 
the low-resolution spectra provided by our SED it is not 
possible to estimate the redshift since only one possible 
emission line is visible. However, the spectrum is clearly 
very blue and the extend of the object is only slightly 
larger than the point spread function. Therefore the can- 
didate #130 could be consistent with a quasar/AGN. For 
a true identification we have to wait for a better spectrum, 
though. The spectral fits to the RXJ0943.0+4701 X-ray 
emission (see previous section) do not provide strong con- 
straints, but these results are certainly not in contradic- 
tion with RXJ0943.0+4701 being an AGN. So combin- 
ing all these different lines of evidence, our conclusion is 
that the strongly variable X-ray source RXJ0943. 0+4701 
is most likely an AGN. 

2.7. The quasar at z=2.055 

In order to determine the X-ray brightness of the quasar 
at redshift z=2.055 (Dressier et al. 1993) we use the same 
method as for RXJ0943. 0+4701: a smoothed model of the 
point spread function (in this case the point spread func- 
tion of the HRI) is subtracted from the image in Fig. 1. We 
find that the quasar has about 30 counts above the clus- 
ter emission in the HRI observation. This corresponds to 
a count rate of 7x 10~ 4 counts/s. Assuming a Galactic hy- 
drogen column density of 1.27x 10 20 cm~ 2 (Dickey & Lock- 
man 1990) and a power law spectrum with photon index = 
2.3, we find a flux /(0.2keV) = 2.7 x 10~ 31 crg/cm 2 /s/Hz 
and a luminosity of log[L(0.2kcV)] = 28.1 (correspond- 
ing to £x(0.1-2.4keV)= 1.4 x 10 45 erg/s), which places 
the quasar among the most X-ray luminous quasars (La 
Franca et al. 1995; Green et al. 1995). However, since the 
quasar it seen through a galaxy cluster, it is unavoidable 
that the gravitational lens effect of the cluster magnifies 
the flux of the quasar. Therefore its true X-ray luminosity 
is probably somewhat lower. 

A comparison of the HRI and the PSPC observation 
(Fig. 6a and b) might suggest a possible X-ray variabil- 
ity of the quasar, because in the PSPC observation Ml is 
slightly brighter than M2 and it is the other way round in 
the HRI image (the latter is smoothed such that the res- 
olutions of both images are equivalent). But the number 
of photons in the two maxima are not significantly differ- 
ent. As these small differences are within the statistical 
fluctuations there is no indication for a brightening of the 
quasar between the two observations. 
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object 


a(2000) 


5(2000) 


R magnitude 


probable identification 


#94 


09 42 59.1 


47 01 04 


20.37 


spiral 


#99 


09 42 58.2 


47 01 02 


19.75 


elliptical/Sa at z=0.30-0.35 


#115 


09 42 55.8 


47 00 56 


19.26 


M0-M2 or K6-K8 star 


#118 


09 42 59.3 


47 00 56 


21.05 


Sb, possibly cluster member 


#130 


09 42 57.2 


47 00 50 


20.60 


blue compact object, slightly larger than PSF 


#134 


09 42 59.0 


47 00 48 


21.63 


possibly elliptical and cluster member 


#158 


09 42 57.5 


47 00 35 


21.61 


Sb-Sc, possibly cluster member 


#166 


09 42 58.1 


47 00 29 


18.48 


M0-M2 or K6-K8 star 



Table 3. List of possible optical counterparts of the X-ray source RXJ0943. 0+4701 and probable identifications. 
Numbers are as in Fig. 8. 



3. Discussion and Conclusions 

The main results from the X-ray analysis of CL 0939+4713 
are summarized in Table 4. It is quite surprising for such 
a rich cluster to have an X-ray luminosity typical for a 
less rich cluster; e.g. its X-ray luminosity is smaller than 
the mean luminosity of the richness class 4+5 clusters in 
Soltan & Henry (1983) and CL 0939+4713 is probably 
richer than these clusters (Dressier et al. 1994). Of course 
there is a large range of luminosities within each richness 
class. Burg et al. (1994) find a factor of 25 between the 
fifth and 95 th percentiles, although there seems to be a 
tendency that the range is becoming smaller when going 
to higher richness classes. As Burg et al. did the analysis 
only up to richness class 2, it is hard to quantify how 
unusual the X-ray luminosity of CL 0939+4713 is. 

One might try to explain the low X-ray luminosity 
by the non-equilibrium configuration of the cluster, which 
is obvious from the pronounced substructure visible in 
the ROSAT/HRI image. Evidently, there are still sepa- 
rate units in this cluster which are still in the process 
of merging. After the final merging the cluster will have 
probably a higher central density and - as the luminosity 
is proportional to the square of the density - also a higher 
luminosity in this stage. But as this effect is not expected 
to be large (Schindler & Mullcr 1993) it cannot explain 
the relatively low luminosity of this obviously pre-merger 
cluster. 

The temperature could be determined here without 
the contamination of the variable source RXJ0943. 0+4701 
to T = 7.6+f^ keV. The luminosity-temperature rela- 
tion for high redshift clusters (e.g. Mushotzky & Scharf 
1997; Tsuru ct al. 1997) predicts a temperature around 6 
keV, i.e. the temperature of CL 0939+4713 is somewhat 
higher than expected from its luminosity, but still con- 
sistent within the scatter of the luminosity-temperature 
relation. From the point of view of the velocity disper- 
sion the temperature is consistent within the scatter of 
the velocity-temperature relation (Edge & Stewart 1991; 



countrate HRI [counts/s] 
Lx(0.1-2.4keV) [erg/s] 
Lx(bol) [erg/s] 
temperature [keV] 
metallicity [solar] 
n H [cm -2 ] 


0.023 
6.4± ( ^ x 10 44 
1.6+°^ x 10 45 
7.6± 2 ; 8 6 
0.22t^ 4 
0.9+g-S x 10 20 


subcluster 


Ml 


M2 


M gaa (r < 265 kpc)[M Q ] 
M tot (r < 265 kpc)[M ] 
gas mass fraction 


4.1 x 10 12 
6.0l 2 'g x 10 13 
7 + 2% 


4.2 x 10 12 
7.4+H x 10 13 
6 + 2% 



Table 4. Summary of the X-ray properties of CL 
0939+4713, and its two subclusters Ml and M2. 



Lubin & Bahcall 1993; Bird et al. 1995), but again slightly 
on the high temperature side. Taking into account the sub- 
structure in the cluster one can explain such a behaviour 
by a merger roughly perpendicular to the linc-of-sight. 
The increase of temperature (and of luminosity) is visible 
from all directions while the increase in the velocity dis- 
persion is very small when it is observed from a direction 
perpendicular to the collision axis. 

It seems puzzling that a cluster with such a high 
galaxy content has not a high metallicity. The metallicity 
of 0.22^Q 22 is in the normal range compared with other 
clusters with this temperature (and less galaxies) which 
have Fe abundances around 0.3 (Fukazawa et al. 1996, 
Mushotzky & Loewenstein 1997). The ratio of iron mass 
to light makes the difference even more evident. The value 
found in CL 0939+4713 of M Fe /L B w 5 x 10~ 3 M Q /L Q 
is smaller than the one expected for rich clusters of 0.01- 
0.02 (Rcnzini et al. 1993), although one must keep in mind 
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Fig. 8. Optical image in the region around RXJ0943. 0+4701. The image has a size of about 80 arcsec by 80 arcsec 
(east is left) and is taken in the I-band filter. The possible optical counterparts of RXJ0943. 0+4701 are marked with 
numbers (see Table 3). The position of the X-ray source RXJ0943. 0+4701 (M3) of the PSPC observation is marked 
with (x) for the soft band maximum, (+) for the hard band maximum and (*) for the broad band maximum. The 
corresponding marginal maximum in the HRI observation is marked with (o). 



the errors are so large that there is some overlap and, in 
addition, the observed dispersion in metallicity is about a 
factor of two (Mushotzky & Loewenstein 1997). 

In order to find an explanation for this difference we 
have a closer look into the galaxy population of the clus- 
ter, having in mind the conclusions of Arnaud et al. (1992) 
that only ellipticals and lcnticulars contribute to the metal 
enrichment of in the intra-cluster medium. Andrcon et 
al. (1997) found from HST images that CL 0939+4713 
is overabundant in spirals. They compared the percent- 
age of ellipticals and lenticulars with the one in the Coma 
cluster. While they find about 50% E's and SO's of identi- 
fied galaxies in CL 0939+4713, 80% of the Coma cluster 
galaxies are E's and SO's. Furthermore, we can compare it 
with another cluster which is not so well virialised as the 
Coma cluster but has also pronounced substructure like 
CL 0939+4713: the Virgo cluster. This cluster has a per- 
centage of ellipticals and lenticulars of 81% (Binggeli et 
al. 1985) - also much larger than in CL 0939+4713. These 
are all results which are based on a morphological classi- 



fication. Using a spectral analysis one finds a somewhat 
different result for CL 0939+4713 (Belloni et al. 1995): 
ellipticals 67%, Spirals and Im 15%, E+A galaxies 19%. 
The high fraction of E+A galaxies would be classified in 
the morphological analysis as ellipticals. But as they have 
a preceding starburst phase, they are in terms of the iron 
enrichment more similar to spirals. Therefore one should 
not count the E+A galaxies as ellipticals, so that the el- 
liptical fraction is just 67%. This number is larger than 
the one from the morphological study but still smaller 
than other elliptical percentages. Summarizing, the low 
elliptical content together with the low iron mass to light 
ratio confirms the correlation of elliptical and iron content 
found by Arnaud et al. (1992). 

The gas mass fraction in this cluster of 6 or 7% is 
lower than in other clusters, where 10-30% are observed 
(Bohringer 1995). The small values can be either due to 
an error in the determination in the total mass, as the 
total mass for non-relaxed clusters is general difficult to 
determine (Ota et al. 1998; Allen 1997). A major prob- 
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mass fraction a compact structure of the subclusters, 
which is found in the subclusters around Ml and M2 
(see small core radii in Table 1). For a standard-shaped 
isothermal cluster the temperature has a fixed relation 
with the total mass, as well as the X-ray luminosity with 
the gas mass. To get a low gas mass fraction a more com- 
pact gas distribution than the average one is necessary, 
which produces the same X-ray luminosity with a smaller 
amount of gas. As this is the case in CL 0939+4713 every- 
thing gives a consistent picture. The mass to light ratio 
of M/Lb ~ 35OM /L is also rougly in the standard 
range. This implies that there is not necessarily an error 
in the mass estimate and the low gas mass fraction can 
be real. Such variations in the baryonic fractions would 
have implications for cosmology. They would be explain- 
able either by dynamical processes operating differently on 
baryonic and dark matter or by pre-inflationary baryon- 
to-total mass fluctuations which are preserve d into the 
post-inflationary epoch (see Evrard 1997). 

This new X-ray analysis of ROSAT/HRI and ASCA 
data confirms that CL 0939+4713 is a young cluster. It 
is a cluster where mergers on different levels are taking 
place, which is obvious from the complex subcluster struc- 
ture. Also other properties point in this direction: for its 
richness the cluster has a relatively low luminosity, low 
gas mass, low total mass and low iron mass, which are all 
indicators for young systems. 

The variable source RXJ0943. 0+4701 at a distance of 
about 1.5 arcmin from the cluster centre was discovered by 
comparing the ROSAT/PSPC and the ROSAT/HRI im- 
age. Its luminosity decreased by at least a factor ten within 
the five years between the two ROSAT observations. Only 
a small fraction of X-ray sources has such large variations 
in the flux (Voges & Boiler 1997). 

To identify RXJ0943.0+4701 we use the X-ray proper- 
ties (X-ray spectrum, short term time variability) in com- 
bination with optical properties taken from observations 
with various filters (SED, ratio of optical to X-ray flux). 
This comparison rules out stars and cataclysmic variables 
as origin of the X-ray emission of RXJ0943. 0+4701, so 
that the most likely origin would be an AGN. 

So far only very few AGNs have been found with 
large X-ray variabilities and these strongly variable AGNs 
turned out to be mainly narrow-line Seyfert 1 galaxies (e.g. 
Piro et al. 1988; Grupe et al. 1995, Otani et al. 1996, Boiler 
et al. 1997). It is not clear yet what physical processes are 
responsible for the strong X-ray variability. Possible expla- 
nations (as discussed in Boiler et al. 1997) could be strong 
relativistic effects and partial covering by occulting struc- 
tures in an accretion disc. Optical follow-up observations 
will give us detailed information on this interesting object. 
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Fig. 3. Top: elliptical f3 model (see Table 1); the contours 
are logarithmically spaced with a spacing of 0.03. The 
central surface brightness is 1.3 x 10~ 2 cts/arcmin 2 /sec. 
Bottom: Residual image after subtracting the elliptical 
(3 model. The image is smoothed with a Gaussian filter 
of a =10 arcsec. The contours represent the significance 
of the emission above the (3 model in units of a. Dot- 
ted contours represent negative values. The positions of 
Ml (square), M2 (triangle) and the quasar at redshift 
two (plus sign) are indicated. The substructure in the 
west around Ml has a significance of more than 4a. It 
is clearly extended and structured. The emission of the 
z=2.055 quasar has a significance of more than 3<r. 
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Fig. 4. X-ray spectra of CL 0939+4713 obtained with the 
ASCA/GIS, ASCA/SIS as well as the ROSAT/PSPC with 
the rcdshift as a free fit parameter (a) and with fixed red- 
shift (b). Each panel shows the GIS spectrum (top), the 
SIS spectrum for the energy range of 1-10 keV (middle) 
and the PSPC spectrum (bottom). The normalisation of 
the PSPC spectrum is shifted arbitrarily for the image. A 
Raymond-Smith model (folded with the response of the 
corresponding instrument) which gives the best fit for the 
GIS, SIS and PSPC spectra simultaneously is superposed. 
The contamination due to RXJ0943. 0+4701 is subtracted 
from the PSPC spectrum. 
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Fig. 7. Lightcurve of RXJ0943.0+4701 from the 
ROSAT/PSPC observation. It is consistent with a 
constant brightness over the observing interval of 72 
hours. The observation was carried out in eight intervals 
with exposure times between 1300 sec and 2400 sec. The 
lightcurve is background subtracted, the cluster emission 
at the position of RXJ0943.0+4701 is not excluded. As 
the cluster emission should be constant it should give 
only an offset to the curve. Because of the limited number 
of photons the error bars are quite large. 



